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Executive  Summary 


The  classic  orthogonal  expansion  of  antenna  patterns  by  spherical  har¬ 
monics  is  finding  new  applications  in  computational  electromagnetics,  an¬ 
tenna  design,  and  wireless  systems.  For  example,  recent  research  uses  spher¬ 
ical  harmonics  to  decouple  the  multiple  antennas  in  a  multiple-input  multiple 
output  (MIMO)  radio  system  from  the  propagation  medium.  This  approach 
reveals  “that  channel  capacity  is  dependent  on  the  number  and  distribution 
of  the  spherical  modes.”  Consequently,  the  spherical  harmonics  force  an 
end-to-end  approach  to  MIMO  systems  designs — adaptively  excite  selected 
spherical  modes  to  maximize  channel  capacity. 

This  report  reviews  spherical  harmonic  expansions  to  compute  wideband 
antenna  patterns.  Some  “antenna”  patterns,  such  as  the  far-field  of  a  ship, 
require  lots  of  computer  time  per  each  frequency.  Consequently,  computing 
such  a  pattern  across  a  wide  frequency  band  at  each  frequency  is  prohibitive. 
Spherical  harmonics  offer  an  alternative — compute  the  pattern  over  a  coarse 
sampling  across  the  frequency  band,  interpolate  the  spherical  harmonics,  and 
use  these  harmonics  to  recover  the  pattern  at  a  fine  frequency  sampling  with 
small  CPU.  This  report  develops  an  adaptive  spline  interpolation  algorithm 
with  error  control.  The  development  used  antennas  with  smooth  analytic 
patterns  for  verification.  Further  developments  must  test  this  algorithm  on 
a  suite  of  tactical  antenna  patterns. 
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1  Spherical  Harmonics  and  Antenna  Patterns 

Spherical  harmonic  expansions  can  speed  antenna  and  electromagnetic  (EM) 
computations.  The  spherical  harmonics  {Y™}  are  an  orthonormal  basis  for 
the  square-integrable  functions  on  the  sphere  [1,  page  572]:  if  g  G  L2(§), 

oo  n 

s(M)  =  £  y  CC (»,*)■ 

n— 0  m=—n 

One  natural  application  for  spherical  harmonic  expansions  is  in  the  ray¬ 
tracing  codes  that  require  the  near-held  antenna  patterns  to  be  integrated 
with  the  ray  tracings.  Figure  1  shows  a  logarithmic  antenna  and  its  patterns 
at  3  and  6  GHz  [2], 


(a)  (b)  (c) 


Figure  1:  (a)  logarithmic  antenna  (b)  pattern  at  3  GHz  (c)  pattern  at  6  GHz 
[2,  Figure  4]. 

If  this  antenna  operates  in  the  indoor  environment  of  Figure  2(a)  and  the 
wideband  channel  response  in  Figure  2(b)  is  to  be  computed,  each  ray  that 
penetrates  the  transmitter  and  receiver’s  patterns  must  be  weighted  with  the 
pattern’s  amplitude  and  phase  over  a  dense  sampling  of  the  frequency  band. 

The  computation  of  all  the  patterns  over  a  wide  and  densely  sampled 
frequency  band  is  considerable  and  the  standard  (/,  9,  0)-table  lookup  slows 
the  ray-trace  computation.  The  spherical  harmonic  expansion  takes  the  form 
[10] 

oo  n 

9(».*/)  =  E  E  ccwry.y 

n= 0  m=—n 

Each  coefficient  g™(f)  is  a  function  of  frequency  that  is  decoupled  from  the 
angles.  More  practically,  the  compact  antenna  need  only  a  few  spherical 
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Figure  2:  (a)  indoor  rays  (b)  time-domain  channel  response;  [2,  Figure  1]. 


coefficients  to  approximate  an  antenna  pattern  [2].  Because  the  spherical 
harmonics  Y™(6,  (p)  admit  fast  computations,  spherical  expansion  avoids  the 
slower  table-lookups.  More  importantly,  the  coefficients  g™(f)  are  smooth 
functions  of  frequency.  If  an  antenna  pattern  requires  a  day  to  compute  per 
frequency,  computing  a  dense  sampling  of  the  patterns  across  the  frequency 
band  is  impossible,  ffowever,  if  the  patterns  can  be  interpolated  with  con¬ 
trolled  error  using  only  a  few  samples,  any  sampling  density  of  the  patterns 
across  the  frequency  band  is  possible. 

This  report  offers  a  simple  adaptive  spline  interpolation  (ASI)  algorithm 
to  benchmark  more  sophisticated  pattern  approximation  methods.  Section  2 
sets  out  two  antenna  patterns  to  be  tested  in  the  spline  algorithm.  Section  3 
reviews  the  spherical  harmonic  functions  Y™(9,  (p)  and  their  expansions.  Sec¬ 
tion  4  uses  one  of  the  antenna  patterns  to  develop  the  ASI  algorithm  and 
the  error  bounds.  Section  5  tests  the  ASI  algorithm  on  the  more  difficult  an¬ 
tenna  pattern.  Section  6  discusses  applications  of  spherical  harmonics  with 
an  emphasis  on  MIMO  wireless  that  include  polarizations. 
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2  Short-Wire  Antennas 


Two  short-wire  antenna  are  reviewed  to  set  the  notation  and  motivate  the 
spherical  harmonic  approach.  Figure  3  illustrates  spherical  coordinates  around 
a  short  antenna  carrying  a  z-directed  current  in  free  space: 

i 


Figure  3:  Spherical  coordinates. 

The  magnetic  vector  potential  is  [4,  Eq.  13-47] : 

0 
0 

z 

When  the  wire  carries  a  uniform  current  Iq  amps  at  frequency  /, 

J {z'\t)  =  z IQ  exp(+ j2vr/t) , 


A(r;  t)  = 


47T 


rh/2  j (z’\t 

' -h/2 


r  (z 


,-i ' 


r \z' 


-dz'- 


r  z  = 


(1) 
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the  magnetic  vector  potential  is 


A(r;  t)  = 


4:71 


'  —h/2 


r-ru 


»+ml,  fk/2  exp(-j2T/c-1||r-rM||})J-, 


=  Z/o4tC 


dz' 


'  -h/ 2 


r-rr 


=  zIn—e+j27vft  [ W 2  exP(-Jfcllr-r(^)ll})^/ 

°4vr  ./_fc/2  ||  r  —  r(V)|| 

where  k  denotes  the  wavenumber  with  wavelength  A: 

27 if  27T 


k  = 


A 


With  a  slight  abuse  of  notation,  the  Fourier  component  of  the  magnetic  vector 
potential — in  cylindrical  coordinates — admits  the  representation 


rh/2  exp  (  -j  k ,  r 2yn  +  (z  -  z' } 2 ) 

A(rcyn,  z,  (f>]  f)  =  z/0-p-  /  -  V  = - J-dz'.  (2) 

J  —h/2  7^2  _l_  ( -y  —  2 


rcyn  +  -  * 


The  held  is  rotationally  constant  around  the  z-axis  and  symmetric  with  re¬ 
spect  to  z.  Figures  4,  5,  6,  and  7  illustrate  the  magnetic  vector  potential  in 
range  and  elevation  as  a  function  of  increasing  wire  length.  For  this  short- 
wire  antenna,  these  plots  show  that  the  number  of  lobes  in  the  pattern  is 
approximately  the  ratio  of  the  wire’s  length  h  to  the  wavelength  A. 
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Magnetic  Vector  Potential  ||  A  ||  :  Short  Wire 


Figure  4:  r-z  slice  of  the  magnetic  vector  potential  (h  =  A). 


Magnetic  Vector  Potential  ||  A  ||  :  Short  Wire 


Figure  5:  r-z  slice  of  the  magnetic  vector  potential  (h  =  2A). 
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Magnetic  Vector  Potential  ||  A  ||  :  Short  Wire 


r  (m):  /  =1  (A),  k=3  (m);  ^=300  MHz 


Figure  6:  r-z  slice  of  the  magnetic  vector  potential  (h  =  3A). 


Magnetic  Vector  Potential  ||  A  ||  :  Short  Wire 


r  (m):  /Q=l  (A),  h=4  (m);  ^=300  MHz 


Figure  7:  r-z  slice  of  the  magnetic  vector  potential  (h  =  4A). 


6 


When  the  current  has  the  more  physic  ally-plausible  distribution  [7,  Eq. 
2-124] 

J  (z';k)  =  zImsm(k(h/2  —  I^D), 

the  magnetic  vector  potential  in  the  radiation  zone  (r  3>  h )  admits  the 
approximation  [7,  page  82] 


A(r ,  9,  &  k)  «  z  Ime 


-jkr  k  cos(khcos(9))  —  cos(kh/2) 


2rrr 


ksin(0y 


(3) 


The  held  is  rotationally  constant  around  the  ^-axis  and  symmetric  with 
respect  to  z.  Figures  8,  9,  10,  and  11  illustrate  the  magnetic  vector  potential 
in  range  and  elevation  as  a  function  of  increasing  wire  length. 
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V  (m):  7^-1  (A);  h~\  (m);  f-300  MHz 


Figure  8:  x-z  slice  of  the  magnetic  vector  potential  with  sinusoidal  current 
( h  =  1  m). 


Magnetic  Vector  Potential  ||  A  || :  Short  Wire 


Figure  9:  x-z  slice  of  the  magnetic  vector  potential  with  sinusoidal  current 
(h  =  3  m). 
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v  (m):  7^-1  (A);  h~ 5  (m);  f-300  MHz 


Figure  10:  x-z  slice  of  the  magnetic  vector  potential  with  sinusoidal  current 
(h  =  5  m). 


Magnetic  Vector  Potential  ||  A  |  :  Short  Wire 


v  (m):  7^-1  (A);  hrl  (m);  f-300  MHz 


Figure  11:  x-z  slice  of  the  magnetic  vector  potential  with  sinusoidal  current 
(h  —  7  m). 
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3  Spherical  Harmonic  Expansions 

The  spherical  harmonics  are  computed  for  a  short-wire  antenna  to  get  a 
baseline  example  for  the  spline  interpolants.  The  spherical  harmonics  {Y™} 
are  the  functions  defined  on  the  sphere  [1,  Eq  12.149] 

U”(«,  0)  =  |c“(cos(#))e^,  (4) 

where  the  P™’ s  are  the  associated  Legendre  functions.  Figure  12  plots  the 
normalized  associated  Legendre  functions  [1,  Eq  12.148] 


Pnm(cos(0))  =  4 
for  n  =  4  and  m  =  0,  1,  2,  3,  4. 


1 2n  +  1  (n  —  m) ! 


(n  +  m)\ 


^(COSW) 


Figure  12:  r-z  slice  of  the  magnetic  vector  potential. 


Although  P™’s  are  orthonormal  with  respect  to  the  inner  product  for 
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fixed  m  [1,  Eq.  12.107] 


(T'^V%)=  /  7T(cosWW(cos((9))  sin (9)dB  =  5n 


it  is  the  complex  sinusoid  exp (jmd>)  that  forces  orthogonality  over  both  in¬ 
dices  [1,  Eq.  12.150]: 


y™(M)^(M)sin(0)d0# 


Cm.  m'  a 


n,n'  • 


The  spherical  harmonics  {Yff}  are  an  orthonormal  basis  for  the  sqnare- 
integrable  functions  on  the  sphere  [1,  page  572]:  if  g  e  L2(E> ) 


oo  n 

9(e,4>)  =  J2  E  iccio,*)- 

n= 0  m——n 


Example  1  (Sinusoidal  Current)  Let  g(f)  €  L2(§)  be  the  collection  of 
patterns 


cos(/c/rcos(0))  —  cos(kh/2 ) 
k  sin(6))2 


(27 xf  =  kc) 


extracted  from  Equation  3  that  are  indexed  by  frequency  f .  All  the  spherical 
harmonics  vanish  for  m^O  because  this  pattern  is  constant  in  (f: 


n= 0 

Because  this  pattern  symmetric  in  z,  All  the  spherical  harmonics  vanish  for 
n  odd  because  this  pattern  symmetric  in  z: 


g(9A-,{)  =  E  9l(/)U°„(M). 

n= 0 

Figure  13  plots  estimate  of  these  coefficients  for  the  short-wire  antenna  with 
1-meter  length  and  the  frequencies  fo  =  300  x  {1,  3,  5,  7}  MHz  (The  lines  are 
added  for  clarity).  The  plot  shows  that  the  coefficients  roll  off  rapidly  with 
increasing  index.  The  plot  also  shows  that  the  higher  frequencies  take  longer 
to  roll  off  in  concordance  with  the  lobing  seen  in  Figures  8,  9,  10,  and  11. 

Figure  If  plots  the  first  few  Legendre  coefficients  as  functions  of  fre¬ 
quency.  These  coefficients  plot  like  the  Bessel  functions  with  the  “DC”  coef¬ 
ficient  g[)(/)  decreasing  as  the  frequency  increases. 
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10 


Short-Wire  .Antenna:  1  meters  vertical  length 


n 

Figure  13:  Legendre  coefficients  for  selected  frequencies. 


MHz 

Figure  14:  Selected  Legendre  coefficients  as  functions  of  frequency. 
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Example  2  offers  two  guides  for  sparse  interpolation: 

•  If  increasing  frequency  always  forces  the  most  detail  in  the  pattern, 
expanding  at  the  highest  frequency  shows  the  number  of  spherical  co¬ 
efficients  to  retain. 

•  If  the  spherical  coefficients  always  show  this  smooth  sinusoidal  function 
of  frequency,  spline  interpolation  will  be  sufficient  for  interpolation. 

The  next  section  develops  an  adaptive  spline  interpolation  using  this  exam¬ 
ple. 
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4  Adaptive  Spline  Interpolation  (ASI) 

The  problem  of  adaptive  interpolation  for  pattern  synthesis  assumes: 

•  computing  a  pattern  g(0,  0;  /)  is  costly, 

•  the  spherical  coefficients  {g™(f)}  vary  smoothly  as  functions  of  fre¬ 
quency, 

•  the  number  of  spherical  coefficients  to  retain  has  been  determined  from 
the  highest  frequency. 

Under  these  assumptions,  the  spherical  coefficients  {g™}  can  be  approxi¬ 
mated  by  high-accuracy  splines  using  relatively  few  frequency  samples.  This 
section  shows  that  a  sequence  of  splines  can  rapidly  approximate  the  spher¬ 
ical  coefficients  with  computable  error  bounds.  A  spline  algorithm 

is  developed  using  the  short-wire  antenna  of  Example  1.  The  algorithm  is 
based  on  the  rapid  convergence  of  the  splines  to  a  smooth  function. 

Theorem  1  (Spline  Convergence)  [11,  Theorem  2. 4.3. 3] 

Let  g  G  CA ( [/min ,  /max])-  Let  A  denote  a  partition 

A  {/min  =  fo  ^  ^  f Np  =  /max} 

with  mesh  size  and  deviation  from  uniformity  measured  as 
IIAIU  =  max{ | fn  -  fn_ i|  :  n  =  1, . .  .,NF}, 

IIAH-oo  =  min{ | fn  -  fn_ i|  :  nf  =  1, . . . ,  NF}. 

If  9a  is  the  spline  that  interpolates  g  on  A  and  satisfies  g'^(f)  =  g\f )  at 
f  =  train ,  /max,  then  there  exists  constants  Ck  <  2  that  do  not  depend  on  the 
partition  such  that 

119a1 -»WIU  <  Ct||S(4»|u]A£f,  k  =  0, 1,2,3. 
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If  the  partitions  are  uniform,  the  splines  converge  as 


ll9i-9lU<C„||9<4>|U||A||i,. 


If  gpjp  denotes  the  spline  computed  from  the  uniform  partition  of  Np  +  1 
frequency  points,  the  splines  converge  as  0[Np 4]: 

ll9»-9lU  <  Qll9|4)ltoo  (/m“jVF/m1")  '  =:  (5) 

That  is,  doubling  the  number  of  sample  points  reduces  the  error  by  a  factor  of 
16  and  reuses  the  previously  computed  sample  points.  This  spline  algorithm 
is  developed  using  the  pattern  g(9,  0;  /)  obtained  from  the  short-wire  antenna 
of  Example  1.  This  pattern  has  the  expansion 


/) 


cos(kh  cos(d))  —  cos(/ch/2) 
k  sin(d)2 


n= 0 


Figures  15,  16,  and  17  show  the  splines  computed  from  5,  9,  and  17  equally 
spaced  frequency  samples  for  a  1-meter  wire  (h  =  1).  The  spherical  coef¬ 
ficients  {.92n (■/')}  are  the  solid  lines.  The  dashed  lines  are  the  splines.  The 
circles  mark  the  interpolation  points  that  are  the  only  information  available 
to  the  analyst.  The  splines  show  rapid  convergence  to  the  unknown  spherical 
coefficients  {g^nif  )}-  Theorem  1  allows  error  estimates  to  be  extracted  from 
this  sequence  of  splines. 
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Figure  15:  Legendre  coefficients  and  five-point  splines. 


MHz:  HI  g,-g  1^=0.066159 

Figure  16:  Legendre  coefficients  and  nine-point  splines. 
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Figure  17:  Legendre  coefficients  and  17-point  splines. 

There  is  the  possibility  of  estimating  the  error  \\g—gNF\\oo  even  though  the 
function  g  is  unknown.  Table  1  compares  these  unknown  errors  \\g  —  g^p  ||oo 
to  the  computable  estimates  \\gNF  —g2NF\\oo-  The  latter  is  a  strong  predictor 
of  the  unknown  error.  The  triple  bar  is  the  sup-norm  error  over  all  the 
coefficients: 

Ills  -  9nf\ Hoc  :=  rnax{||  g™  -  g™Nf\\oo  :  n  =  0, 1, . . . ,  N;  -n  <  m  <  n}. 

The  remainder  of  this  section  develops  two  error  estimates  and  incorporates 
these  estimate  into  the  adaptive  spline  algorithm. 

The  first  error  estimate  is  obtained  from  the  observation  that  “doubling 
the  sample  points  reduces  the  error  by  a  factor  of  16.”  This  approximation 

||s  S2iVj?||oo  ~  txIIs  Sw^Hoo- 
16 

and  the  triangle  inequality 

IIS  —  SiVF||oo  <  ||S  —  92Nf\\oo  +  \\Q2Nf  —  9Nf\\oo 

~  TX  1 1 S  9nf  1 1  oo  T  ||  92NF  9nf  1 1  oo 


0 

So 

0 

0 

*4 

0  - 

£(, 

V — Gk — — 3 

v-r  y 
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Table  1:  Spline  convergence. 


NF 

III  9  fitypHloo 

|||<?2iVF  —  9nf  oo 

15  1 1 92Np  QNj?  |  oo 

CF 

4 

0.5251 

— 

0.5683947 

— 

8 

0.066159 

0.53287 

0.0696491 

129 

16 

0.0028367 

0.065296 

0.0029790 

252 

32 

0.00012164 

0.0027928 

0.0001219 

173 

64 

0.0000098425 

0.00011428 

— 
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yield  the  “16/15”  estimate  of  the  error  listed  in  Table  1: 

16, 


F||oc  ^  *7/VF||oo* 


\\9~9n. 

The  second  error  estimate  uses  Inequality  5 


\\g  —  9nf  ||oo  < 


CF 

AT4 


to  estimate  Cf  from  the  spline  over  the  finer  mesh: 

CF  CF 

II  9Nf  -  <72JVf||oo  <  II  9Nf  -  9 1 1  oo  +  II 9  -  ^JV^Hoo  <  + 


N%  '  161V4' 

This  inequality  produces  a  lower  bound  for  CF  listed  in  Table  1: 

/p 

■|| 9nf  ~  92Nf ||oo  <  Cf. 


16  Ni 


17 


(6) 


(7) 
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Figure  18  applies  these  estimates  to  the  data  in  Table  1.  The  unknown 
error  1 1  <7  —  gNp  \ \  x  is  plotted  as  the  blue  line  marked  with  circles.  The  “16/15” 
estimator  of  Equation  6  is  plotted  as  the  dashed  red  line.  This  estimator 
closely  tracks  this  error  but  cannot  predict  the  error  for  arbitrary  Np.  The 
Cp/Np  estimator  is  plotted  as  the  dotted  green  line  using  Cf  =  252  obtained 
from  Np  =  16.  This  estimate  predicts  that  Np  =  64  is  necessary  to  push  the 
error  below  ICE4.  As  the  number  frequency  points  double,  the  analyst  can 
monitor  the  running  error  and  update  the  Cp/Np  estimate. 
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This  process  of  doubling  the  number  of  frequency  samples  and  monitoring 
the  convergence  of  the  splines  is  summarized  in  the  following  algorithm.  At 
each  stage,  the  running  error  estimate  is  computed  by  the  “16/15”  estimate 
Ci.  The  predicted  number  of  sample  points  is  computed  from  the  asymptotic 
bound  62  and  reported  for  CPU  decisions.  When  the  running  estimate  beats 
the  user’s  error  bound,  the  last  computed  spline  is  reported  and  its  error 
bound  is  estimated  using  the  asymptotic  bound  e2. 

Algorithm  1  (Frequency  Doubling) 


1.  Estimate  the  number  of  Legendre  coefficients  to  compute. 

2.  Set  the  error  e. 

3.  Set  NF  =  4  and  compute  . 

(a)  compute  g2NF- 

(b)  estimate  a (NF)  :  =  (16/15)  x  \\g2NF  -&vF|l°o- 

(c)  compute  Cnf  from  Equation  7. 

(d)  estimate  e2(NF)  max{ G'nf } /Nf4  and  reportthe  estimated  sam¬ 
pling  necessary  for  e2(NF)  <  e. 

(e)  break  if  €i(Nf)  <  e  or  set  NF  :=  2 NF  arid  loop  back  to  (a). 

4.  Report  the  approximation  as  the  last  spline  gNF • 

5.  Estimate  the  error  as  e2(W)- 
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5  Testing  the  ASI  Algorithm 

When  the  wire  carries  a  uniform  current  Jo  amps  at  frequency  /, 

J  (z'\t)  =  z/oexp(+j2vr/f), 

the  magnetic  vector  potential  is  computed  by  Equation  2  and  plotted  in 
Figures  4,  5,  6,  and  7.  The  magnetic  vector  potential  in  the  radiation  zone 
(r  h )  at  the  frequency  of  the  current  [7,  Eq.  2-122]  is 

pjkr  rh/2 

A^r.M;/)  «  —  /  J(z'yk‘  '“><»>  dz' 

4vr r  J  ~h/2 

jkr  rh/2 

=  I0—  ejkz  cos(0)  dz' 

4vr r  J_h/ 2 

e jkr  rh/2 

=  Iq -  /  cos(kz'  cos(#))  dz' . 

4tt r  J_h/2 

For  comparison  with  the  sinusoidal  current,  define  the  uniform  current  pat¬ 
tern  as 

rl/2 

gu{9,4>',  f)  :=  /  cos{hz'  cos(0))  dz' . 

J -1/2 

This  section  applies  the  Frequency-Doubling  Algorithm  to  estimate  the  spher¬ 
ical  coefficients  {g™}  of  the  expansion 

OO  Tl 

=  £  V  gZ(f)YZ(f>AY 

n= 0  m=—n 

Example  2  (Uniform  Current)  Because  this  pattern  is  constant  in  <f>  and 
symmetric  in  z, 

OO 

n=0 

Figure  19  plots  numerical  approximations  of  these  coefficients  for  a  1 -meter 
length  and  the  frequencies  f0  =  300  x  (1,3, 5,  7}  MHz  (The  lines  are  added 
for  clarity).  The  plot  shows  that  the  coefficients  roll  off  rapidly  with  increas¬ 
ing  index.  Figure  If  plots  the  first  few  Legendre  coefficients  as  functions  of 
frequency. 
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Figure  19:  Legendre  coefficients  for  selected  frequencies. 
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Figure  20:  Uniform  current  (selected  Legendre  coefficients  as  functions  of 
frequency). 
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Figure  21:  Uniform  current  (Legendre  coefficients  and  five-point  splines). 
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Figure  22:  Uniform  current  (Legendre  coefficients  and  9-point  splines). 
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Figure  23:  Uniform  current  (Legendre  coefficients  and  seven-point  splines). 
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Table  2  reports  on  the  convergence  of  Algorithm  1.  The  triple  bar  is  the 
sup-norm  error  over  all  the  coefficients: 


III  9  ~  9Np  1 1 1  oo  :=  max{||5f™  -  g™Nf\\O0  :  n  =  0, 1, . . . ,  TV;  -n  <  m  <  n}. 

This  error  would  not  be  known  to  the  analyst  because  the  coefficients  are 
not  known.  This  error  can  be  using  the  “16/ 15” -rule  estimates  the  error  of 
spline  gNp  using  the  adjacent  spline  (]2Nf ■  Table  2  reports  that  this  estimate 
(column  4)  is  bounding  the  error  in  (column  2)  For  example,  the  error  in 
spline  g16  is  estimated  using  g32  as  the  proxy  for  g: 


III#  —  016  | 


16 

15 


1 1  #32  —  #16  | 


0.0005352. 


If  the  analyst  terminates  with  NF  =  32,  Algorithm  1  estimates  the  error  in 
#32  the  coefficients  CF  in  column  5: 


III#  #32 1 


cv 

naf 


—  =  0.00003147. 
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Table  2:  Spline  convergence  for  the  uniform  current. 


nf 

III#  9Np  1 1 1  oo 

1 1  |#27Vp  —  gNF  oo 

jf  ||#27Vf  -  #VF|  oo 

cF 

4 

0.1343 

— 

0.1440 

— 

8 

0.004521 

0.1350 

0.004763 

33 

16 

0.0005190 

0.004465 

0.0005352 

17 

32 

0.00003824 

0.0005017 

— 

30 
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6  Applications 

This  section  reviews  recent  applications  of  spherical  harmonic  expansions 
to  multi-antenna  systems,  antenna  bounds,  and  exploiting  the  sparsity  of 
the  spherical  pattern  expansions  with  compressed  sensing.  The  spherical 
harmonics  separate  the  interaction  of  the  antenna  and  the  communication 
channel  and  shows  that  the  harmonics  can  be  interpreted  as  “modes”  to 
connect  the  transmitter  to  the  receiver  [5],  CY:2010,  [6].  These  topics  extend 
the  capabilities  of  antenna  performance,  design,  and  wireless  propagation  at 
SPAWAR. 


6.1  MIMO 


This  report  stated  with  the  applications  of  the  spherical  harmonic  expansions 
to  indoor  wireless  using  a  single  transmitter  and  receiver.  Figure  24  is  a 
system  diagram  of  a  MIMO  radio  using  two  antennas  at  the  transmitter  and 
two  antennas  at  the  receiver. 


MIMO  Channel 


Figure  24:  2x2  MIMO  system  schematic. 


The  received  signal  s n(t)  is  modeled  as  the  convolution  of  the  transmitted 
signal  Sx{t)  with  the  channel  matrix: 


The  channel  matrix 


OO 


H{t)  =  Y,Hn(t)6{t-Tn), 


n=  1 


26 


includes  all  the  rays,  fading  processes,  and  antenna  gains  in  each  random  pro¬ 
cess  {Hn(t)}.  Glazunov’s  2009  Pli.D  thesis  made  a  fundamental  application 
of  spherical  harmonics  to  MIMO  systems  where  he  shows  that  [5]: 

. . .  optimum  decorrelation  of  the  antenna  signals  is  obtained  by 
the  excitation  of  orthogonal  spherical  vector  waves  modes. 

That  is,  optimal  MIMO  performance  can  be  computed  using  these  modes  [6]: 

.  . .  the  expansion  of  the  MIMO  channel  matrix,  H,  in  spherical 
vector  wave  modes  of  the  electromagnetic  held  of  the  antennas  as 
well  as  the  propagation  channel . . .  [that]  maps  the  modes  excited 
by  the  transmitting  antenna  to  the  modes  exciting  the  receive 
antennas  and  vice  versa. 

This  approach  computes  MIMO  capacity  [13]: 

. . .  the  use  of  spherical  harmonics  . . .  decouples  the  effects  of  the 
antennas  and  the  propagation  medium.  . . .  [and]  illustrate  how 
channel  capacity  is  dependent  on  the  number  and  distribution  of 
the  antennas  spherical  modes 

This  approach  computes  the  capacity  of  a  MIMO  system  that  uses  polar¬ 
ization  as  part  of  the  MIMO  diversity.  Figure  25  is  a  schematic  of  a  2  x  2 
MIMO  system. 


Figure  25:  2x2  dual-mode  MIMO  [9]. 

From  [9] :  “The  four-arm  antenna  structures  provide  orthogonal-polarized 
operational  modes  by  feeding  one  of  the  two  adjacent  antenna  arm  pairs  with 
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an  odd  phase  relation.”  The  spherical  mode  expansion  can  accurately  model 
the  antenna  correlations.  Because  antenna  correlation  is  directly  related  to 
the  MIMO  channel  capacity  or  diversity  gain,  the  spherical-mode  expansions 
provide  a  modal  approach  to  analyze  the  MIMO  channel. 

6.2  Antenna  Gain 

The  classic  antenna  bounds  of  Chi  have  been  generalized  by  Kogan  [8] : 

The  algorithm  of  calculation  of  a  high  bound  of  antenna  gain 
occupying  area  of  space,  restricted  some  sphere  and  operating  in 
given  frequency  band  is  given.  It  is  supposed,  that  the  antenna 
is  linear,  obeys  to  a  principle  of  a  causality  and  passivity.  The 
algorithm  is  based  on  determination  of  a  extremal  gain  of  an 
antenna  unit  radiating  a  sole  spherical  harmonic. 

Kogan’s  approach  synthesize  Chi’s  circuits  with  the  wideband  Fano  compu¬ 
tations  using  the  spherical  harmonics  to  link  a  circuit  and  a  mode.  This 
approach  may  apply  to  bound  the  limits  of  wideband  patterns. 

6.3  Compressed  Sensing  and  Pattern  Interpolation 

If  a  signal  is  “sparse”  in  the  frequency  domain,  then  sub-Nyquist  sampling 
in  the  time  domain  can  recover  to  the  signal  to  arbitrary  accuracy.  Likewise, 
for  those  antenna  patterns  that  have  sparse  spectral  harmonics,  the  antenna 
patterns  can  be  recovered  from  a  limited  number  of  samples  on  the  sphere 
[12]. 

This  compressed  sensing  result  will  not  reduce  the  method  of  moment 
computations.  All  the  current  must  be  computed  before  any  part  of  the  an¬ 
tenna’s  pattern  can  be  calculated.  However,  an  application  may  be  found 
in  an  inverse  problem:  recovering  and  identifying  a  pattern  from  a  limited 
number  of  observations.  Consequently,  there  may  be  a  natural  connection  be¬ 
tween  general  theory  of  compressed  sensing  and  the  specific  radon-transform 
algorithms  of  ISAR. 
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